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In this paper we analyze the throughput and delay characteristics 
of the Carrier-Sense Multiple-Access protocol with a queueing model. 
The effects of finite buffer size, bursty arrivals, and collision detection 
with exponential rescheduling are examined. 

The Carrier-Sense Multiple- Access protocol could be used on a bus 
in a packet switch. It works by sending a signal to all sources when 
the bus is occupied. A source postpones transmission to the bus when 
this signal is heard. Since the signal takes a positive amount of time 
to reach the sources, two or more sources occasionally will attempt to 
use the bus at about the same time. When this occurs, all the packets 
are destroyed and rescheduled. 

We conclude that: (i) it is important to choose the mean resched- 
uling time correctly, and (ii) performance degrades significantly 
when compound Poisson arrivals (peaked traffic) replace Poisson 
arrivals (smooth traffic). 



I. INTRODUCTION 

The Carrier-Sense Multiple-Access (csma) protocol has been pro- 
posed for resolving conflicts when several sources attempt to use a 
single channel. In this paper we investigate the throughput and delay 
characteristics of csma. 

1.1 Background 

When several sources attempt to use a single channel, a protocol (in 
this context, protocol is synonymous with queue discipline) is required 
to allocate the channel among the sources. When the channel is 
occupied, a signal is sent indicating that state of affairs. When a source 
wants to use the channel, it first listens for this signal. If the signal is 
not heard, the source starts transmitting. If the signal is heard, the 
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source postpones transmitting and tries again at another time. The 
advantage of this system is that a device to control the sources is not 
required. The disadvantage is that occasionally messages will be de- 
stroyed because two sources will transmit at the same time. This is a 
consequence of the fact that signals travel at finite speeds, so there is 
a delay between the epoch when one source seizes an idle channel and 
the epochs when the other sources can first hear the busy-channel 
signal. Thus, soon after a source seizes the channel, another source 
may sense that the channel is free even though the channel is busy. 
When this occurs, both messages are destroyed because their bits have 
been merged. At the end (earlier if collision detection is used) of each 
transmission, the source determines whether the transmission was 
successful. If it was, the source goes about its business; if it was not, 
each message is rescheduled as if the channel were sensed as busy. 

The csma protocol is used in the Ethernet* local-area distribution 
system and has been considered for part of a digital switch. Since csma 
is envisioned as a protocol for packet networks, we will refer to the 
arrivals as packets and the channel as a bus. 

1.2 Relation to other work 

The first study of csma is Kleinrock and Tobagi. 1 In addition to the 
version of csma studied here, which they call nonpersistent csma, they 
investigated various forms of persistent csma and slotted csma that 
are not alluded to in this paper. Their paper implicitly assumes that 
the rescheduling delay is infinite. This means that packets which find 
the bus occupied or which are destroyed in a collision are abandoned. 
For the most part, we assume, as their paper does, that packets arrive 
according to a Poisson process. In a subsequent paper by Kleinrock 
and Tobagi, 2 the rescheduling times are modelled as geometrically 
distributed random variables. They used a finite number of sources 
and assumed the arrivals to be quasi-random (i.e., finite-source Pois- 
son). Also, they introduced an unnatural discretization of the time 
scale, which causes some small distortions in the results. A continuous 
version of this model is described and solved in Halfin. 3 Our analysis 
shares many features with the analysis found in Ref. 2, particularly 
the exploitation of the regeneration epochs. 

Tobagi and Hunt 4 add collision detection to the model of Kleinrock 
and Tobagi. 2 Collision detection is a feature that informs a source that 
its packet has been destroyed soon after the collision occurs. Our 
model of collision detection is based on the model presented in Ref. 4. 

Rappaport 5 and Rappaport and Bose 6 describe an elaborate model 
of a protocol similar to csma. 



* Ethernet is a trademark of Xerox Corporation. 
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1.3 Summary of results 

The models considered in this paper can be easily solved numeri- 
cally. The following conclusions were reached by considering many 
numerical examples. 

(i) It is important to choose the mean rescheduling time correctly; 
performance is not sensitive to changes in this number near its best 
value. 

(ii) If the mean rescheduling time is based on a designed load but 
the realized load is different, the realized performance is almost as 
good as if the mean rescheduling time were based on the realized load. 

(Hi) When the traffic intensity is no larger than 0.7, the throughput 
is insignificantly lower than the traffic intensity. 

(iv) Performance is significantly degraded when compound Poisson 
arrivals (peaked traffic) replace Poisson arrivals (smooth traffic). 

(v) When the traffic intensity is no larger than 0.7, collision detec- 
tion lowers the average waiting time significantly and reduces the 
sensitivity of the throughput to the mean rescheduling time. For higher 
traffic intensities, collision detection significantly increases the 
throughput and decreases the average delay. 

1.4 Outline of this paper 

Our model is described in detail in Section II, and the method of 
solution is outlined. The details of the solution are presented in Section 
III. Numerical examples and empirical conclusions from the model are 
given in Section IV. In Section V, we consider the effects of bursty 
traffic by introducing compound Poisson arrivals. We return to Poisson 
arrivals in Section VI and examine the benefits of collision detection. 
Appendix A records the transition probabilities omitted from the text, 
and Appendix B lists the most important symbols used in the text. 

II. MODEL AND OVERVIEW OF THE SOLUTION METHOD 

Our setting is a queue of the M/D/l/K type. Let X be the rate of the 
Poisson arrival process. There is a single server with a constant service 
time per packet. At most K packets may be in the system (queue plus 
server) at any time. Packets that arrive when the system is full are lost 
forever and have no effect on the system. Our use of constant service 
times reflects the assumption that all packets have the same length. 

The parameter K has two potential interpretations. One is to suppose 
there is a buffer that can hold K - 1 packets, so the maximum number 
of packets in the queue is K - 1. The other is to suppose there are K 
ports and a packet that cannot seize a port is lost. We will show that 
proper limiting probabilities will not exist without this bound on the 
number of packets in the system. 
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When a packet is on the bus (i.e., when it enters service), a signal is 
sent that warns other packets that the bus is occupied. (Equivalently, 
a signal indicating the bus is free is turned off when a packet is on the 
bus.) Let h denote the one-way propagation delay for this signal. It is 
common to all potential users. The propagation time between a source 
and the bus depends on the distance between them. If we calculate h 
for the source that is farthest from the bus, the assumption that h, the 
one-way propagation delay, is common to all users will underestimate 
the performance measures. This propagation delay has two effects. 
The first effect is that for the first h time units after a packet occupies 
a previously idle bus, other packets that arrive and want to use the 
bus will do so because the signal that the bus is occupied has not 
reached them. When this happens, we say that a "collision" has 
occurred and the first h time units of a service time are called the 
"vulnerable" period. The second effect is that the signal that the bus 
is free is not received until the bus has been free for h time units. 

When a collision occurs, the bits in both packets are scrambled, and 
so both packets must be retransmitted. We assume that the collision 
is detected at the end of the unsuccessful transmission. Then each 
source involved in the collision draws a number from an exponential 
distribution with mean 1/a, which is the length of time until the source 
next attempts to use the bus. Furthermore, all sources that attempt to 
transmit a packet and hear the signal that the bus is occupied make 
their next attempt in this fashion. We call a the "retry rate." 

Our objective is to express the throughput (equivalently, the asymp- 
totic departure rate or the asymptotic proportion of time that the bus 
is successfully transmitting packets) and average waiting time of a 
packet (i.e., the departure time minus the arrival time) to the system 
parameters A, K, h, and a. 

2. 1 Outline of solution 

Choose time units so that the processing time of a packet is unity. 
Consider a packet, 0, that gets on an empty bus at time t. This pac- 
ket will relinquish the server at t + 1 and the signal that is on the 
bus will cease being sensed at some random time t with t + l + AS 
t < t + 1 + 2h. If no packets arrive during the vulnerable period, then 
r = t+l + 2h. 

Suppose we assume that 2P holds the bus for a (nonrandom) length 
of time v, with 1 < v < 1 + 2h, at the end of which ^relinquishes the 
server and the signal that 0* is on the bus ceases. Retries that occur 
during the vulnerable period destroy the messages involved in the 
collision. These retries are not considered to have seized the bus. 

When v = 1, these assumptions will provide an upper bound for the 
throughput and a lower bound for the average number of customers 
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present (in a buffer or on the bus) in the steady state. When v — 1 + 
2h, a lower bound for the throughput and an upper bound for the 
average number of customers present is obtained. Little's theorem (the 
queueing formula L = XW) implies that setting v = 1 provides an 
upper bound for the average waiting time and that setting v = 1 + 2h 
provides a lower bound. Therefore, we will solve a model with constant 
service times p. 

Since h is typically much smaller than one (we use h = 0.01 in our 
examples), we expect that the bounds will be close together. That is 
the case in our numerical examples. Therefore, either of the bounds is 
a good approximation to the true measure of performance. 

Let X(t) be the number of packets present at time t. We want to 
obtain 

Pi = \imP{X(t) = i) (1) 



t~ 



for each i = 1, 2, • • • , K. These "limiting probabilities" are easily 
shown to be independent of the distribution of the initial state X(0). 
To obtain {/>,}, we embed a discrete Markov chain at customer 
"ejection" epochs, which are defined as follows. Let C m be the epoch 
where the mth packet to get on an empty bus leaves the bus. This may 
be the end of a successful transmission where the packet leaves the 
system or it may be the end of a destroyed transmission where the 
packet rejoins the queue. To encompass both types of events, we call 
C m the mth ejection epoch. 

Let Y m = X(C m ) A lhn si o X(C m + s); it represents the number of 
customers present just after C m . Since the arrivals are Poisson and 
retrials are governed by an exponential distribution, it is easy to see 
that { Y m \ m = 0, 1, • • • } is a Markov chain. It is also easy to see that 
this chain is irreducible and aperiodic (details are omitted) so that the 
limits 

•m A lim P{Y m = i), i - 0, 1, • • • , K, 

m— »oo 

exist and are independent of Yo, are positive, and sum to one. 

Three steps are used to obtain {pi} from {777}. The first step is to 
relate {tt,} to {pj} where 

p\ A lim P{Y n = i and transmission is successful}, 

71— »00 

i-0,1, ...,JC-1. (2) 

The relationship is obtained by calculating the state-dependent prob- 
ability of having a collision. The second step is to use a "rate up equals 
rate down" argument to show 

Apffpi, i = 0,l,... ,K-l, (3) 
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where f is the ejection rate. The third step is to obtain p by a renewal- 
reward argument. Then eq. (3) is used to calculate f and the remaining 
Pi's. 

III. SOLVING THE MODEL 

In this section we give details of the solution of the model described 
in Section II. 

3. 1 Collision probabilities 

Let c„ be the probability that n packets arrive in an interval of 
length, v\ then 

c n = e~ Xv — — , n = 0, 1, ... . 
n\ 

A well-known property of Poisson arrival processes (e.g., Corollary 5- 
13 in Heyman and Sobel 7 ) is that the arrival epochs are iid and uniform 
over [0, v] when it is given that n > arrivals occurred during [0, v\. 
Thus, for a service interval that starts with i packets in the queue, we 
have 

Vn A P{no arrivals prior to h\n arrivals in a 

service interval that starts with i in queue} 

n = 0, 1, • • • and any i. (4) 



Use the memoryless property of the exponential distribution to obtain 

5, A P{no retries prior to h\n arrivals in a 
service that starts with i in queue} 

= e- M , i = 0, 1, • • • and any n. (5) 

Since a collision is avoided if, and only if, there are no arrivals and no 
retries during the vulnerable period, we have 

d„(i) = P{no collision and n arrivals in a service 
interval | start with i in queue} 



= VnCnSi, i, n > 0, (6) 



and 



d n {i) A ^{collision and n arrivals in a service 
interval | start with i in queue} 

- (1 - Vn8i)c n , h n > 0. (7) 

Equations (6) and (7) are used to calculate the transition probabilities 
of the embedded Markov chain { Y m ; m = 0, 1, • • • } . 
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3.2 Transition probabilities 

Recall from Section 2.1 that Y m is the number of packets in the 
system just after the mth ejection epoch. All of these packets must be 
in the queue because of the propagation delay in broadcasting that the 
bus is available. Let 

Pu±P{Y m+i =j\Y m = i). (8) 

In this section we present formulas for computing py, < i,j < K. It 
will be convenient to call both exogenous packet arrivals and retries 
"arrivals"; the former are "outside" arrivals. Then, for 2 < i < K — 1 
and < i + n < K - 2, we may write 

Pu+n = P{ next arrival is from outside | i in queue} 

X [P{n outside arrivals during service and no collision | service 
starts with i in queue} 

+ P{n — 1 outside arrivals during service and a collision | service 
starts with i in queue}] 

+ P{next arrival is a retry | i in queue} 

X [P{n + 1 outside arrivals during service and no colli- 
sion | service starts with i — 1 in queue} 

+ P{n outside arrivals during service and a collision | service 
starts with i — 1 in queue}]. 
Use the memoryless property of the exponential distribution to obtain 

Pfnext arrival is from outside \i in queue} =- — . 

A + ice 

Set /? = A/a for notational simplicity; then 

Pi, i+n = t-^ [dnd) + dn-i(i)] + t^-z [^„ + i(i - 1) + d n (i - 1)]. (9) 

To obtain the remaining entries of (p,y) we need to account for 
boundary conditions. The details and the formulas are given in Ap- 
pendix A. 

It is easy to see that (/?,>•) is irreducible and aperiodic. Therefore, 
there is a unique stationary distribution that is also the limiting 
distribution. Furthermore, since the continuous- time process {X{t)\ 
t > 0} starts from scratch (i.e., it regenerates) each time the Y-process 
reaches zero, the X-process is regenerative and the mean time between 
regeneration points is finite. This implies the limits in eq. (1) exist. A 
similar argument establishes that the limits in eq. (2) exist. 

We can use eq. (9) to explain why the finite capacity, K, is essential 
for our analysis. Suppose K = co; then eq. (9) is valid for i > 2 and 
every n > 0. Equation (9) is also valid for n = — 1 if we set dj{i) = 
whenever / < 0. Define 
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gi = X npi,i+ n — pa-it 

1 

it is the expected size of a jump out of state i in the Y-process. 
Intuitively, if linv-,. gi > 0, the Y-process is tending to infinity and 
lim m _coP{Y m <y'} =0 for ally. From eq. (9) we can calculate (details 
are omitted) 

/?(1 - e- iah e-^) i[e- a{i - l) (e- Xv - e~ xh ) + e^e - ^ 

8i = Xv+ TTp TTp ; 

and so, if h > 0, 

lim gi = Xp. (10) 

This phenomenon can be described physically. Wheny* is large, e~ Jah 
is small so the probability that an ejection is a completion is very 
small. This means that once state j is reached (as it must be because 
the process is irreducible), the number of customers present grows 
monotonically with very high probability. 

The corollary in Kaplan 8 states that an irreducible and aperiodic 
Markov chain for which eq. (10) holds is not ergodic if p,> = whenever 
j < i — k for some k that is independent of i. At most, one departure 
can occur at a time, so p„ = whenever j < i - 1. Therefore, the Y- 
process is not ergodic when K = oo. 

3.3 Solving the balance equations 

Instead of solving ir = irP to obtain {in}, we will solve balance 
equations between two sets of states. Let 

K 

Pi(>:j) 4 P{ Y m+ i >y| Y m - i} - S Pik. 

In the steady state, the rate of transitions between states {0, 1, • • • , i} 
and states {i + 1, i + 2, • • • , K} must be the same in both directions, 
so* 

i 

m+iPi+u = S iTkPk{>i + D, i = 0, 1, • . • , K - 1. (11) 

We can solve eq. (11) by replacing 777 by Xi and setting x = 1, and using 
eq. (11) to obtain Xt+i from x , • • • , Xi. By setting 

l K 

m = Xi X X J 
I y-o 

we obtain a nonnegative solution of eq. (11) that sums to one. 



* The basic idea was developed by Robert Morris and Eric Wolman. The details for 
Markov chains are given in Theorem 7-13 of Heyman and Sobel. 7 
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3.4 Departure point probabilities 

Recall that C m is the mth ejection epoch. Let us write C m = D when 
C m is a departure epoch. Now 

P{C m +i = D,Y m+i =j) 

= 2 P{C m+ i = D, Y m+1 =j\ Y m = i)P{Y m - 

for every / and m. We will see below that P{C m +\ = D, Y m+1 = 
j\Y m = i) does not depend on m, so denote it by qy. Then 

K 

p'j = lim P{C m +\ = D, Y m+ \ =j] = 2 qi/ni, 

m-.oo ,=0 

7 = 0,1, ...,#-1. (12) 

Observe that, for 1 <j < /iT - 2 and 1 < i <j - 1, P{C m+1 = D, Y m+i 
= J\Y m = i) 

= P {next arrival is from outside | i in queue} 

X P {j — i outside arrivals during service and no collision | service 

starts with i in queue} 
+ P {next arrival is a retry | i in queue} 
X P {j — 1 + 1 outside arrivals during service and no collision | service 

starts with i — 1 in queue} 

- dj-i(i) + -^—. dj-i + i(i - 1). 



/? + i J (S + i 

The remaining entries are obtained by accounting for boundary be- 
havior. The results are: 

^oo = Cotjo = Co, (13a) 

"-TT7T (13b) 

and 

q i0 = for i>2. (13c) 

Forl<y<ii:-2, 

<7o, = dj(0), (14a) 

qtj = -^—. dj-id) + — l — dj- i+1 (i - 1), 1 < i <>, (14b) 
P + i P + i 

qj+1J = /+j+l <? ° (,/) ' (14C) 
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and 

g«,--O f iay + 2. (14d) 

When 7* ■■ K - 1 we obtain 

Gur-i - S <*„(0), (15a) 

P + 1 tf-l-i P + * X-i 

l<i<^T-l (15b) 

Qkjc-i " 8/c-i- 

In eq. (46) we show that the infinite sums in eq. (15) have represen- 
tations as finite sums containing no more than K terms. Since C m +\ = 
D and Y m+i = K cannot occur simultaneously, que ~ for every i. 
Notice that (the subscript oo denotes a limit) 

nek i p;-P{c.-Ay w sjr-i} 



= P{an ejection epoch is a departure epoch} 

= P{no collision}, (16) 

where the last two probabilities are steady-state quantities and can be 
interpreted as long-run proportions. 

3.5 Converting {p!} to {p ( } 

Let E(t) be the number of ejections by time t. Since {E(t)\ t > 0} 
regenerates when an ejection leaves the system empty (i.e., when 
Y m -= for some m), 

f A Urn E(t)/t 

exists. It is the ejection rate. It is also the rate at which packets gain 
access to an idle bus, so f is the arrival rate of new and rescheduled 
packets, which is denoted by G in Kleinrock and Tobagi. 1 

The rate at which the number of packets present jumps from i to 
i + 1 is Ap,. The rate at which it jumps from i + 1 to i is £pj. In the 
steady state these rates must be equal so we have 

Lemma 1: Ap, - ft»J, i - 0, 1, - • . , K - 1. (17) 

Lemma 1 can be proved rigorously. The proof uses standard methods 
and is omitted. 
There are K + 1 unknowns in eq. (17) and K equations. We will find 
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p by an independent argument; the remaining pfa and f are obtained 
from eq. (17). 

Let To be the amount of time X( ■ ) is zero in an arbitrary cycle of the 
X-process, and let M be the expected length of a regeneration cycle. A 
basic property of regenerative processes (see, e.g., Theorem 6-7 in 
Heyman and Sobel 7 ) is that 

po - E(To)/M. (18) 

Since X(t) =0 if, and only if, t is in an idle period, 

E(To) = 1/A. (19) 

When Y = i, let \pi be the average time between en ejection epoch and 
the next ejection epoch. We have 



*- 



v + 



ia + X 



if i<K 



„ + if i = K, 

Ka 



(20) 



because »//, is v plus the expected time to the next arrival after the 
ejection epoch. Let m, be the mean number of visits of the Y-process 
to state i during an arbitrary regeneration cycle of the X-process. Thus, 
m, is the mean number of visits of the Y-process to state i between 
visits of the Y-process to state zero, because the X-process regenerates 
whenever the Y-process enters state zero. It is easy to show (see 
Exercise 7-78 in Heyman and Sobel 7 ) that 



mi = 777/770, i = 0, 1, 
Using eqs. (20) and (21) we obtain 

v , I: (,■ + -' 

,=0 no 



* 1 

M = J turn = 



K-\ 

1 1Ti[v + -. 

\ la + X 



,K. 



+ "* ,,, + E 



(21) 



1 1 



K-\ 



TTo A V + 0+i 



fai (5-JTK 



+ 



K 



(22) 



where p = Xv. 

Substituting eqs. (19) and (22) into (18) yields 

po = wo / ( p + l im + ~z~ 



(23) 



P + i ' K 
Obviously [and formally from eqs. (12) and (13)], we find that 

Po = 7To. (24) 
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Combining eq. (17) with i = and eqs. (23) and (24) yields 

Equation (25) shows that f can be computed when the balance equa- 
tions are solved. 

3.6 The throughput, occupancy, and average waiting time 

The throughput, 0, is the asymptotic departure rate. Since the 
departure process regenerates whenever the ejection process regener- 
ates, 6 is well defined. The asymptotic departure rate equals the 
asymptotic rate at which packets are accepted, because all accepted 
packets eventually depart and the number of packets present is no 
larger than K\ therefore 

= A £ Pi- (26) 

o 

Combining eq. (26) with eqs. (16) and (17) yields 

6 = fa. (27) 

This equation states that the departure rate equals the ejection rate 
multiplied by the asymptotic proportion of ejections that do not suffer 
a collision. Equation (27) illuminates the essential trade-offs involved 
in using csma. One expects that £ increases and n c decreases as A and 
a increase. This means that for each A, there is a value of a that 
maximizes 6. Let 0* (A) be the largest value of that can be achieved 
when A is specified. We also expect that 0*(X) will first increase with 
X and then start decreasing. 

Let <j> be the long-run proportion of time that the bus is occupied. 
The regenerative arguments used above can be used to prove that <f> is 
well defined and that 

K 

$ - v £ m-i/M. (28) 

o 

Using eqs. (22) through (25) in conjunction with eq. (28) yields 

<t> - H. (29) 

This equation can be obtained from Little's theorem by regarding 
the bus as "the system." Then £ is the arrival rate (since only packets 
that depart could have arrived) and v is the expected time a packet is 
on the bus. Little's theorem asserts that \v is the average number of 
packets on the bus, which is the proportion of time that the bus is 
occupied. Combining eqs. (27) and (29) yields 

<{t=v0/n c . (30) 
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Suppose that we attempt to increase 6 by increasing the arrival rate, 
and K and a are adjusted to keep n c constant (i.e., a fixed proportion 
of collisions). Equation (27) shows that increases in 6 must be accom- 
panied by decreases in the occupancy of the bus, and that doubling 
the throughput would halve the occupancy. 

From eqs. (17) and (25), p, for 1 < i < K - 1 is obtained in the 
obvious way and p K is obtained from p K = 1 - Yfl~* P>- Then 

K 

L = £ ipi 

i 

is the average number of packets present in the steady state. Since 6 
is the arrival rate of packets that enter the system, Little's theorem 
yields 

W=L/6, (31) 

where W is the average length of time that a packet is in the system. 

3.7 Relation to the model of Kleinrock and Tobagi 

In this section we explain how the "basic equation for the through- 
put" [eq. (3) in Ref. 1] can be obtained from the model in this paper. 
This will clarify the similarities and differences between the two 
models. 

In Ref. 1 Assumption 1 states that the average retry interval is large 
compared with the packet transmission time. In our notation, this 
assumption is that 1/a is large compared with p. Since there is no 
parameter corresponding to a in eq. (3) of Ref. 1, it appears that they 
have set a = 0. Let us do that. This means that every transmission 
that is destroyed by a collision stays in the queue forever. Conse- 
quently, Assumption 2 in Ref. 1, which states that the interarrival 
times of the point process consisting of packet arrival epochs and retry 
epochs is a Poisson process, is valid because there are no retry epochs 
and the arrival epochs form a Poisson process. 

Another consequence of a = is that without a finite buffer, 
lim/_ooP{X(0 = i) = for all i because each packet transmission has 
a positive probability of being destroyed by a collision. Therefore, we 
make these two assumptions. 

Assumption 3. There is a finite buffer that can hold K — 1 > 
packets. 

Assumption 4. If a packet transmission is destroyed when the buffer 
is full, that packet is flushed from the system. 

The next assumption is required to replicate the collision process in 

Ref. 1. 

Assumption 5. Packets that arrive when the buffer is full and a 
packet is being transmitted will destroy that transmission. 
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With these assumptions and a = 0, the arguments and formulas in 
Sections 3.2 through 3.5 can be used to obtain the results given below, 
but we will bypass that roundabout route and give a direct argument 
(which is essentially the argument in Ref. 1). 

In the steady state, whenever a packet seizes the bus, there are 
K — 1 packets in the buffer. The mean time between entries to state 
K — 1 is, via eq. (20), v + 1/A, and 1/A is the mean length of stay in 
state K — 1, so 

P^-^Ta-TTu- < 32 > 

v + 1/A 1 + Xp 

The probability that a transmission is not destroyed by a collision is 
e~ xh , so 

e = Xp K -ie- Xh . (33) 

Equation (7) of Ref. 1 and a = show that 

v= 1 + 2h- (l-e- Kh )/X. (34) 

Substituting eqs. (32) and (34) into eq. (33) yields 

Xe- Xh 



d = 



A(l + 2h) + e 



-\h> 



which is eq. (3) in Ref. 1. 

We can conclude that eq. (3) in Ref. 1 is valid when a = and the 
derivation in Ref. 1 requires assumptions 3, 4, and 5. Our model does 
not use assumptions 4 and 5 so it should produce a greater throughput 
when all other factors are the same; this is demonstrated numerically 
in Section IV. Since K = 1 is allowed, this derivation explains why eq. 
(3) in Ref. 1 becomes the single-server Erlang loss formula when there 
are no collisions, i.e., why h { yields eq. (9) in Ref. 1. 

IV. NUMERICAL RESULTS 

A Fortran program to solve the equations in Section III was imple- 
mented in double-precision arithmetic on a PDP-11 70. Because of 
memory limitations and the design of the program, JiT < 30 was 
required. The running time on a crt display terminal with a 9600-baud 
line is less than a twinkling of an eye. 

In all the numerical examples we use h = 0.01. This is the value 
used by other authors. It is a reasonable value to use for a system such 
as Ethernet with loops of about 300 meters carrying 128-byte packets 
at 10 megabits/sec. In general, 

u lr 
sc 
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where 

I = loop length in meters, 

r = transmission rate in bits/second, 

s = packet size in bits, and 

c = speed of light in m/second. 

4. 1 The upper and lower bounds are close 

In Table I we present various cases that demonstrate that the lower 
and upper bounds for 6 and W are close. The subscript "avg" stands 
for average and the entries are obtained by setting v — 1 + h. In all 
these cases, K = 20. 

Table I suggests the approximations 

Wwg = l MWi b + Wub) and aV g = l M0ib + Bub). 

We make no claim for the general validity of these approximations. 
Based on the encouraging results in Table I, and on several dozen 
other examples not reported here, we will henceforth use the case v = 
1 + h as an approximation; there will be no explicit mention that our 
numerical results are approximate. 

4.2 The effects of changing a 

The first two pairs of columns in Table I show that for fixed X, 
different values of a yield different values of the performance measures. 
In this section we explore the consequences of changing a. 

When K and A are given, the maximum achievable throughput is 
attained when there are no lost service times from collisions and 
instantaneous retries, i.e., by an M/D/l/K queue. Let the throughput 
of the M/D/l/K queue be denoted by max . When K = 20 and X = 0.7, 
max = 0.700 to three decimal places. In Table II we show that good 
choices of a will achieve = 0.699, but the throughput for poor choices 
of a have a much lower value. In particular, very small values of a do 

Table I — Evidence that the bounds for 6 and W are close 



A 
a 


0.7 
0.01 


0.7 
3.0 


0.5 
0.5 


0.9 
1.0 


1.0 
1.0 


3.0 
2.0 


dub 
Oavg 

Bit 


0.459 
0.457 
0.455 


0.673 
0.667 
0.660 


0.818 
0.812 
0.806 


0.803 
0.796 
0.788 


0.798 
0.790 
0.782 


0.560 
0.555 
0.549 


w„ 6 

"avg 


42.1 
41.9 
41.7 


11.4 

10.2 

9.1 


19.6 
19.1 
18.7 


19.9 
19.3 
18.7 


22.7 
22.3 
21.9 


35.8 
35.5 
35.1 


avg n c 
avg (J> 


0.991 
0.466 


0.828 
0.814 


0.923 
0.888 


0.861 
0.933 


0.824 
0.947 


0.570 
0.984 
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not do well. This confirms the observation given in Section 3.7 that 
the model in Ref. 1, in which a = 0, underestimates the throughput. 

The explanation for the qualitative properties shown in Table II is 
as follows. When a is very small, a packet that upon arrival finds 
another packet in service will spend a long time waiting to retry. This 
makes collisions rare but tends to keep the buffer positions full, which 
makes the throughput and bus utilization low. As a is increased, the 
retries are more frequent, and although there are more collisions, this 
is more than offset by the shorter times spent waiting to retry. When 
a is made sufficiently large, the retries occur so rapidly that too many 
packets are destroyed and performance degrades. 

A significant feature of the data in Table II is that 9 does not vary 
much for 0.5 < a < 1.6. The variation within the interval 0.5 < a < 1.4 
is in the fourth decimal place. The value of a that achieves the lowest 
value of W also achieves a high value of 0, but the value of a that 
achieves the highest value of 6 has a value of W about 9 percent larger 
than the best value of W we have found. 

Let 6 * denote the largest value of that is found for fixed values of 
A, h, and K. Numerical results not reported here have shown that 0*/ 
#max increases as A decreases. This is not surprising because lower input 
rates result in fewer collisions. Since 8* = 0.6993 and max = 0.700, 6*/ 
#max is very close to one. This suggests that csma will provide good 
throughput performance when A < 0.7, and that collision detection 
schemes will not increase throughput very much when A < 0.7. 

4.3 The effects of changing A 

Take K = 20 and h = 0.01. Choose a so that the throughput is 
maximized. The effects of changing A are shown in Table III. 

Notice that 6 increases with A for A < 2 and then decreases very 
slightly when A = 3. This indicates that the phenomenon of lower 
throughput with higher arrival rate will not occur in the normal 
operating range of A < 1. Table III also shows that a should not 
increase as A increases. This property also appeared in every example 
we tried. 



Table II 


1 — Performance 


measures vs. 


a when K 


= 20 and A 


= 0.7 


a 


0.001 


0.01 


0.1 


0.5 


0.8 


1.0 


e 

w 

Dc 


0.362 
53.9 
0.993 

0.368 


0.457 
41.9 
0.991 
0.468 


0.660 
22.8 
0.979 
0.681 


0.6989 
8.34 
0.968 
0.729 


0.6993* 
6.51 
0.963 
0.734 


0.6992 
6.51 
0.963 
0.734 


a 


1.4 


1.6 


2.0 


3.0 


4.0 


5.0 


e 

w 

n,. 


0.6986 
5.53 

0.949 
0.743 


0.6980 
5.52 
0.943 
0.747 


0.696 
5.87 
0.927 
0.758 


0.667 
10.2 
0.828 
0.814 


0.556 
24.1 
0.612 
0.917 


0.423 
42.1 
0.437 
0.977 
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The average waiting times shown in Table III are about the same as 
the corresponding quantities for 1000 sources shown in Fig. 3 of 
Ref. 2. 

Suppose we design a system for a nominal arrival rate and from time 
to time the actual arrival rate differs from the nominal arrival rate. In 
Table II we see that performance suffers if a poor value of a is chosen. 
Now we show that moderate changes in A will not cause much 
degradation in the throughput rate or the average waiting time. 

We assume that K = 10 and h = 0.01 are held fixed. The nominal 
arrival rate is 0.7 and a = 1.6 yields the largest throughput rate. We 
also assume that when A changes, the new value is maintained long 
enough to ensure that steady-state performance measures are ade- 
quate. If we achieve good performance with arrival rates Ai and A2, it 
is reasonable to suppose that we will achieve reasonable performance 
when A is changing from Ai to A2 . 

We will use the notations 6* and W* for the best values of 6 and W 
we have found by varying a. The entries marked 6 and W are for 
a = 1.6. 

Table IV shows that, for the values chosen, deviations from the 
nominal load will not cause serious performance degradation if a is 
kept fixed. 

4.4 The effects of changing K 

As K is increased, more packets can be stored, so the throughput 
should increase. This cannot be carried too far because as K -* 00, 
lim„_oo Y„ = 00 (as shown in Section 3.2), which causes the probability 
of a collision to approach 1 and -* 0. We have not found a value of 
K to demonstrate this numerically. 



Table III — Performance measures vs. A when K = 20 and h = 0.01 

A 07 09 L0 2.0 3.0 



a 0.8 0.6 0.5 0.5 0.4 

0.699 0.813 0.817 0.818 0.817 

ma » 0.700 0.898 0.975 1.00 1.00 

W 6.51 18.9 21.4 23.9 24.1 

n, 0.963 0.911 0.914 0.913 0.905 

4> 0.734 0.901 0.904 0.905 0.912 



Table IV— 6 and W vs. i 


\ when a = 


1.6, K = 


10 


, and h = 


= 0.01 


\ 0.5 


0.6 


0.7 


0.8 




0.9 


1.0 


e 0.500 
e* o.5oo 

W 2.30 
W* 2.18 


0.599 
0.599 
3.07 
2.65 


0.692 
0.692 
4.37 
4.18 


0.764 
0.764 
6.22 
6.09 




0.801 
0.801 
8.06 
8.05 


0.812 
0.815 
9.39 
9.38 
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In Table V below we use the following notations: 
* = largest value of found 
otg = value of a that achieves 0* 
W* = smallest value of W found 
a w — value of a that achieves W* 
W{a* e )= W when a = a* e 

0{<xw) = when a = otw 

#max = throughput of M/D/l/K queue. 

In Table V, 0* and W* increase with K; a w and a £ decrease with K; 
aw^. cte, and the difference is small and gets smaller as K increases; 
W(ag) is not much bigger than W*; and 0(aw) is not much smaller 
than 0*. 

4.5 Conclusions 

From the data presented in this section, and from dozens of unre- 
ported sets of calculations, we conclude that: 

(i) The bounds on and Ware close, and choosing v = 1 + h yields 
a good approximation. 

(ii) It is important to choose a good value of a, and performance is 
not sensitive to changes in a near the best value. 

(Hi) The value of a that minimizes W is close to the value of a that 
maximizes 0. 

(iv) The best value of a decreases as A and K increase. 
(v) For A < 0.7, 0* is essentially the same as max when K > 5. 
(vi) Although a should vary with A, and W do not significantly 
differ from their best values if a is held fixed while A varies over a 
moderately wide interval. 

V. BATCH ARRIVALS 

In this section we replace the assumption that packets arrive ac- 
cording to a Poisson process with the assumption that packets arrive 

Table V — Performance measures vs. K when A = 0.9 and h = 0.01 



K 


5 


10 


15 


20 


30 


0m„ 


0.842 


0.885 


0.895 


0.898 


0.900 


e 


0.771 


0.801 


0.810 


0.813 


0.814 


8(0$) 


0.771 


0.798 


0.808 


0.811 


0.814 


aS 


3.0 


1.4 


0.8 


0.6 


0.4 


0$ 


4.0 


1.8 


1.0 


0.7 


0.4 


W* 


3.61 


8.05 


13.2 


18.8 


30.5 


W(ol) 


3.66 


8.12 


13.4 


18.9 


30.5 



2040 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1 982 



according to a compound Poisson process. In Fuchs and Jackson, 9 
statistical analysis of call arrival times are given. Two of their conclu- 
sions are as follows: The exponential distribution is a reasonably good 
approximation of the time between bursts, and the size of a burst 
(measured in various ways) has a geometric distribution. The purpose 
of this investigation is to find out how sensitive the performance 
measures are to the assumption of Poisson arrivals. We will see that in 
this model, the throughput can be significantly lower with bursty 
arrivals than with Poisson arrivals with same rate. 

Specifically, we assume that bursts arrive according to a Poisson 
process with rate \ h , and the bursts B U B 2 , • • • are iid with 

P{£i - i} - (1 - £)£'. i-1,2, .-.. (35) 

It would be nice to interpret this process as one where messages arrive 
according to a Poisson process with rate Xb and the y'th message 
consists of a random number of packets with a geometric distribution. 
Unfortunately, we cannot do that in the context of the model described 
in Section II. This arrival process does not conform to the finite buffer 
interpretation of the model because only the current lead packet in 
each message would compete for the bus; we, on the other hand, 
assume that all packets in the buffer compete for the bus. This arrival 
process does not conform to the interpretation that there are K ports 
because that would require that if the first packet (i.e., the message) 
finds a free port, then all the packets in the message would enter the 
system. We assume that only those packets that find a port gain access 
to the system. 

5.1 Details of the arrival process 

Let A(t) be the number of arrivals during (0, T). Since A(t) is a 
compound Poisson process, we know that there are constants M and 
V, such that 

E[A(t)] = Mt 

and 

Var[A(*)] = Vt 
with V > M. We call M the arrival rate. From eq. (35) we obtain 

E(B i )= — — and Var(£i) = (1 _g } 2 - 
Standard calculations yield 

E[A(t)]-\btE{Bi)~YZ-t 
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and 

Var[A(t)] = XbtVariBi) + \ b t[E(Bi)] 2 
_ \bt(£+l) 
(1 - £) 2 ' 
In an obvious way we obtain 

2M 2 , V-M 

Letting z = V/M > 1 yields 

2M , M z - 1 
\ b = - and £ = -. (36) 

1 + 2 2+1 

Equation (36) relates the parameters we might obtain from measure- 
ments, M and z, to the parameters of the model, Xb and £. Let 

6(n,7) = P{Bi + B 2 + • • • Bj - n}; 

it is the probability that n packets are contained in y bursts. The sum 
of iid geometric random variables has a negative binomial distribution, 



b(n,j) = 



if n<j 



for all nonnegative integers y and /i. We can compute 6(/i, /) by 
recursion from 

6(0, 0) = 1, 6(1, 1) = 1 - £, 

HjJ) = (1 - Obij - 1,7 - 1), / > 2, 
and 

b(j + k,j)= J+ " ffly + A-lj), A>1 and y>0. 

5.2 Changes in the Po/sson model 

In this section we describe the changes in the equations in Section 
III that are caused by compound Poisson arrivals. 

Let Cb(j) be the probability thaty batches arrive in an interval of 
length v; therefore, 

C6( y ) = e -v^. (37) 

Let 

d n (i,j) = P {no collision and n arrivals in a service interval \j 
batches arrive and start with i in queue} ; 
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then 

d n (i)= £ d n {i,j)c b (j). (38) 

y-0 

Expand on the argument used to obtain eq. (6) to conclude that 

d n (i,j) = SiT}jb(n,j). (39) 

Similarly, 

dn(i) = £ (1 - Si-qjWnJMj). (40) 

j-o 

The transition probabilities described in Section 3.2 and Appendix 
A are given in terms of the d and d functions, so they do not change 
form. The infinite sums, used in Section A.l of Appendix A, have to be 
recomputed; this is done in Section A.2 of Appendix A. 

The steady-state balance equations depend only on the transition 
probabilities, so everything in Section 3.3 is still valid. No changes are 
required for the equations in Section 3.4. 

Since the arrivals do not form a Poisson process, time-average 
probabilities need not equal arrival epoch probabilities and we cannot 
use the arguments in Section 3.5 to obtain the steady-state probability 
that i packets are present. This means we have not obtained L and W 
for this model. We can obtain the throughput by using the arguments 
in Section 3.5. 

Let £ be the ejection rate andp = lim,_oo.P{X(£) = 0}. In the steady 
state, Xtpo is the rate at which transitions leave state and frro is the 
rate at which transitions enter state 0. Since these rates are equal, we 
have 

£ = Xtpo/ira. (41) 

Equations (18) through (25) are valid for compound Poisson arrivals 
when X is replaced by X b . The throughput is obtained from eq. (27). 

5.3 Numerical examples 

Let 0, be the throughput when z is the variance to mean ratio. In 
Tables VI and VII we use the value of a that maximizes 6, and h = 

0.01. 

In Tables VI and VII we see that batch arrivals can significantly 
degrade throughput, and the degradation increases with z and de- 
creases (except, possibly, for very large A) with X. 

In Section IV we saw that 6 is strongly influenced by a. That 
observation suggests that the throughput values in Tables I and II 
might improve if a varied with z. Our numerical experience suggests 
that, in the vicinity of the best a, Q z is not sensitive to changes in a. 
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Table VI — Throughput vs. M for K = 5 and h = 0.01 



M 


0.1 


0.5 


0.7 


0.9 


1.0 


2.0 


5.0 


a 


20.0 


4.0 


4.0 


3.0 


3.0 


2.0 


2.0 


ft 


0.100 


0.497 


0.667 


0.771 


0.798 


0.828 


0.819 


ft 


0.069 


0.382 


0.535 


0.663 


0.709 


0.819 


0.813 


Os 


0.035 


0.205 


0.304 


0.403 


0.451 


0.741 


0.804 


fc/ft 


0.689 


0.769 


0.868 


0.860 


0.888 


0.989 


0.993 


ft/ft 


0.349 


0.412 


0.456 


0.523 


0.565 


0.895 


0.982 




Table VII— 


-Throughput vs. 


M for K = 


10andh = 0.01 




M 


0.1 


0.5 


0.7 


0.9 


1.0 


2.0 


5.0 


a 


10.0 


1.5 


1.6 


1.4 


1.2 


0.9 


0.6 


ft 


0.100 


0.500 


0.692 


0.801 


0.815 


0.821 


0.816 


ft 


0.069 


0.400 


0.590 


0.737 


0.775 


0.816 


0.805 


06 


0.035 


0.235 


0.371 


0.512 


0.576 


0.798 


0.788 


02/0! 


0.690 


0.800 


0.853 


0.920 


0.951 


0.994 


0.987 


05/01 


0.350 


0.470 


0.536 


0.639 


0.707 


0.972 


0.966 



The largest improvement that we could achieve by changing a was 
0.003. 

VI. COLLISION DETECTION 

Suppose that a time units after the vulnerable period ends, the bus 
is examined for a collision. When a collision is detected, the packet 
transmission is aborted and that packet is returned to the buffer. 
Collision detection reduces the time spent transmitting garbage, so it 
will increase throughput and reduce delays. 

6. 1 Changes in the model without collision detection 

Two changes in the equations in Section III are required to describe 
collision detection. The first change is in eq. (7). Let 

rj" = P{no outside arrivals prior to h\n outside arrivals 
prior to a + h)\ 

then 



tfHjJ + ^l' »-0,l, .... (42) 

Let c" be the probability than n outside arrivals occur in an interval of 
length a + h, i.e., 

A-**"**** n = 0,l,...; 

7l! 

then 

d°(i) 4 P{collision and n outside arrivals by a + h\ service 
starts with i in queue} 
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= (l-- n a n 8 i }c a n , *,/i>0, (43) 

where 5, is given by eq. (5). 

The equations in Sections 3.2, 3.3, and 3.4 are given in terms of the 
functions d and d, so no changes are required. In Section 3.5 we need 
to change eq. (20) because the expected length of time that a packet 
is on the bus is not v. Let 

p c (i) = P{a collision occurs | service starts with i in the queue}. 

There are no collisions if, and only if, no arrivals (inside or outside) 
occur within a time interval of length h. Thus, 

l-p c (i) = 8icS 

or 

p c (i) = 1 - 5,cg, i - 0, 1, • • •, K- 2. (44a) 

When i = K — 1, no outside arrivals are permitted, so 

p c (K-l) = l-8 K - 1 . (44b) 

Let Ti be the transmission time of a packet, given that i packets 
were present at the last ejection epoch, and let vi = E(Ti); then 

_ _ f 1 + h if no collision 
' 1 a + h if a collision ' 



and so 

A 



vi = — 



- {(1 + h)[l -pc(i)] + (a + h)p c (i)) 



la + 

+ ^Y {(1 -»- h)[l -p c (i - 1)] + (a + h)p c (i - 1)} 
la + A 

!_•_!, #i . food) + ipctt - 1) . _, _ , 

-l + A-(l-a) r— -r , 1 = 0, 1, ...,J5T- 1. (45) 

We replace the p in eq. (20) by p„ which induces some obvious changes 
in eqs. (22) through (25). 

6.2 Numerical examples 

In Tables II and III we saw that when A = 0.7, the throughput in the 
vicinity of the best a (0.8) is very close to 0.7; but when a < 0.1 or 
a > 3.0, the throughput can be much smaller than 0.7. With collision 
detection, with a = 0.02, throughputs of about 0.7 can be achieved with 
a as large as 5.0. Collision detection reduces the average waiting time 
by about one-third. This example suggests that when a throughput 
close to the maximum possible can be achieved without collision 
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Table VIII — Performance measures vs. A when 
K = 20, h = 0.01 , and a = 0.02 



\ 


0.9 


1.0 


2.0 


3.0 


a 


4.5 


3.0 


2.5 


2.5 


8 


0.891 


0.935 


0.943 


0.942 


W 


8.53 


14.9 


20.6 


20.9 


n,. 


0.718 


0.673 l 


0.628 


0.619 


* 


0.910 


0.958 


0.969 


0.970 



Table IX — Performance measures without 
collision detection 



A 


0.9 


1.0 


2.0 


3.0 


a 


4.5 


3.0 


2.5 


2.5 





0.434 


0.569 


0.606 


0.602 


W 


44.2 


33.3 


32.4 


32.9 


nc 


0.444 


0.585 


0.625 


0.620 


* 


0.989 


0.983 


0.980 


0.981 



detection, collision detection will lower the average waiting time and 
make the throughput less sensitive to a. 

Now we investigate the effects of collision detection when the 
throughput is significantly smaller than the arrival rate. 

By comparing Tables III and VIII we can see the effects of collision 
detection. Throughput increases significantly for each A and W de- 
creases. An indirect effect of collision detection is that larger values of 
a are best. This reduces the time spent waiting for a retry, makes the 
probability of no collision smaller, and yields a larger occupancy for 
the bus. If the values of a shown in Table VIII were used without 
collision detection, performance would degrade significantly, as shown 
in Table IX. From the last two rows of Table IX we deduce that the 
bus wastes a lot of time transmitting packets that have been destroyed. 
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APPENDIX A 
Transition Probabilities 

In Appendix A we record the transition probabilities omitted from 
the text. 

A.1 Transition probabilities for Section 3.2 

Let 

Sri/, i) - £ d n (i) 

n—j 

and 

S d (j,i)= 1 d n (i). 

n-j 

The former is the probability that more than 7' - 1 packets arrive in a 
service interval and there are no collisions when i packets are in the 
buffer at the start of the service interval. The latter is the correspond- 
ing probability when there is a collision. For computational purposes 
it is important to represent these infinite sums as finite sums. Observe 
that 

Sd(j, i) = 5,7(0, i) - 2 d n (i), j > 1, 

n=0 

and 

SAj, i) = SAO, i) - S d n (i), j a 1; 

then 

&<o, <) = i &<••> - « ■ £ (^r) .- { ^f - «*"*■ < 46 > 

n=0 \ " / "• 

Since Sj(0, i) + S d (0, i) = 1 is easily established, 

S d (0, i) - 1 - e"**; (47) 

this equation also can be obtained by calculating the sum explicitly. 
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To modify eq. (9) when i = we have to delete terms with negative 
arguments. When i + n = K - 1 v/e need to recognize that state K - 
1 is reached \f n> K — \ — i packets arrive and there is no collision. 
When i + n = K, we need to recognize that state K is reached if, and 
only if, at least K — 1 — i packets arrive and there is a collision. 

The following transition probabilities are obtained: 

Po.o = d (0) 

Po.j = 3j(0) + rfy-i(O) 1 sy < K - 2 

po.x-1 - 2 5„(0) + cfe- 2 (0) = firtff- i, o) + d K -2(0) 

K-\ 

Po. K = 2 d n (0) = S rf (A--l,0) 

/r-i 

Pi.. = y4o *d) + YVR [Jl(0) + rfo(0)] 

2<7<A"-2 



Aw = y^-g W* " 2, 1) + tfr-ad)] 

[Sj(JSr-l,0) + djr- 2 (0)] 



l + fi 



Puc =YTp Sd{K ~ 2 > 1} + TTb sak ~ 1,oh 

For 2 < i < A" - 2, 

p.,»-i --r-r- ^5o(«- 1) 

Pw = rf/? ^ + TTfl [Ji(l ~ 1) + Mi " 1)] 

Pv = ^p-g Bw(») + rfy-.-i(0] + j^— [dj- i+1 (i - 1) 
+ dj-i(i -1)], i + 1 < / < A - 2 , 
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p tJ c-i = -r^ [S^iK -l-i, i) + d K -2-i(i)] 

+ ^-^ [Sj{K - i, i-l) + dg-i-Ai - D] 

i + p 

Pi .K = J?-; S d (K - i - 1, i) + ^-;S d (K - i, i - 1) 
i + /? i +p 

K — 1 
Pk-\jc-2 = K _ 1 - g do(K - 2) 

'' S rf -(0, tf - 1) + /"] [&(1, ff - 2) 

A. — L -r p 



h ! -S rf (l, iiC-2) 



Pk-\,k-i 


K 


-1 + /? 


Orfvw 






+ d (K - 


-2)] 


Pk-\,k 


K 


-l + P 


S d (0 


Pk.k-x 


= 8 K 


-l 




Pk,k 


= 1 - 


- 8k-\- 





K-1 + /3 



A.2 The effect of batch arrivals 

Equations (46) and (47) have to be modified for the compound 
arrival process. From eqs. (38) and (39) we obtain 

oo oo n 

StiO, i) = £ d„(i) =n d n (ij)c b (j) 

n=0 n=0 y'-O 

= £ C b (j) £ d n (ij) 

= £ c h (j) £ &%&(», » 

7=0 n=j 

-«,£c 6 (y) £ (— ) &(n,y). (48) 

y=0 B-y \ P / 

To evaluate the double sum, let 

n 

a n = £ b(n,j)c b {j); 
j=o 

it is the probability that n customers arrive during a service interval. 
Let 

i( 2 ) = £ z n a n = £ 2" £ b(n,j)c b (j). (49) 

„=0 n=0 y-o 
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Comparing eqs. (48) and (49) we see that the double sum in eq. (48) is 
A(-) evaluated at z = (v — h)/v. 

Standard generating-function arguments yield 



where 



and 



A(z) = C[B(z)], (50) 

B(z) = Y z k P{B 1 = k} = 



. An ^ » *d-0 



*=o 



i-iz 



C(z) = £ z J c b (j) = e- x ^- z) . 

Substitution into eq. (50) yields 

A(z) = exp(\ b p-^A. (51) 

Evaluating eq. (51) at z = (v — h)/v and substituting the result into eq. 
(48) yields 

— Xbvh 



Sd(0, i) = S.exp 
As before, 



v - {v - h)i 

S d (0, i) = l- Sj(0, i) 



is obtained easily. 



APPENDIX B 
List of Symbols 

The following is a list of symbols and their definitions as used in this 
paper. 
a = retry rate 

h = one-way propagation delay 
K = system capacity in packets 
A = Poisson arrival rate 
M = compound Poisson arrival rate 
n c = P{no collision} 
v = service time (constant) 
</> = proportion of time bus is occupied 
= throughput 

#max = throughput of an M/D/l/K queue 
W = average wait (in queue plus in service) of a packet 
z = variance to mean ratio of compound Poisson process 
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f = ejection rate of packets from the bus. 

An asterisk on a symbol means the best value found. The subscripts 
lb, ub, and avg stand for lower bound, upper bound, and average, 
respectively. 
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